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Abstract. We used high quality BeppoSAX data of 6 Seyfert galaxies 
to test realistic thermal Comptonization models. Our main effort was to 
adopt a Comptonization model taking into account the anisotropy of the 
soft photon field. The best fit parameter values of the temperature and 
optical depth of the corona and of the reflection normalization obtained 
fitting this class of models to the data are substantially different from 
those derived fitting the same data with the power law + cut-off model 
commonly used. The two models also provide different trends and cor- 
relation between the physical parameters, which has major consequences 
for the physical interpretation of the data 



1. Introduction 

The broad band X-ray spectra (2-300 keV) of Seyfert 1 galaxies are generally 
well fitted by a cut-off power law continuum + reflection This model depends 
on three parameters: the spectral index T, the high energy cut-off E c and the 
reflection normalization R. The physical interpretation of these parameters is 
generally done in the framework of thermal Comptonization mechanism, which 
is commonly believed to be the origin of the X-ray emission of Seyfert galaxies. 
Using approximate relations (cf. section 2), it is then possible to derive, from 
values of T and E c , values of the temperature kT e and optical depth r of the 
"comptonizing" hot plasma (the so-called corona). 

If such approximations are sufficient in isotropic geometries, where Comptoniza- 
tion process produces roughly cut-off power law spectra, strong discrepancies 
may appear in anisotropic ones, especially for small optical depth and large tem- 
perature. The contribution of the first scattering order to the outgoing flux is 
then significantly reduced producing a spectral break (the so-called anisotropy 
break) in the spectrum. Consequently the physical parameter obtained fitting 
this class of models to broad band X-ray spectra are substantially different from 
those derived fitting the same data with the power law + cut-off model. In a 
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recent paper (Petrucci et al., POO), we have applied accurate Comptonization 
models to fit high quality data of the BeppoSAX long look at NGC 5548 and we 
have underlined the importance of the anisotropic effects. In the same spirit of 
POO, we extend our application of accurate Comptonization models to a larger 
sample of objects. We report in the following the results of this study. 



2. Data and model fitting 

From the complete sample of Seyfert 1 observed by BeppoSAX we have selected 
six objects (NGC 5548, IC 4329A, NGC 4151, ESO 141-G55, Mkn 509 and NGC 
3783) whose observations have high signal-to-noise ratios and hard spectra yield- 
ing a good detection in the PDS instrument. Concerning NGC 4151, we have 
used the observations of July 1996 and January 1999. 

We have fitted the data using two different models for the primary continuum: 

(1) an exponentially cut-off power law plus a reflection component from neu- 
tral material (PEXRAV model of XSPEC, Magdziarz & Zdziardski, 1995) and 

(2) a thermal Comptonization spectrum from a disk+corona configuration in 
slab geometry (code of Haardt (1994, hereafter H94). The fit parameters of 
H94 are the temperature of the corona kT e , its optical depth r, the tempera- 
ture of the disk kT^ (assuming a black body soft emission) and the reflection 
normalization R. On the other hand, the PEXRAV continuum depends only 
on 3 parameters: the e-folding energy of the cut-off power law E c , the photon 
index T and the reflection normalization R. The temperature kT e inferred from 
the PEXRAV fits is simply computed as kT e = E c /2, keeping in mind that 
such approximation roughly holds for r<l. For t>1, kT e = E c /3 would be 
more appropriate. Knowing the temperature, the spectral index derived from 
the PEXRAV fit can be used to determine the optical depth using the following 

1/2 3 

— - (Shapiro, Lightman & Eardley, 
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relation T — 1 ~ 

4 fcT e r(l + r/3) 

1976). This equation is valid for r > 1, and we have checked a posteriori that 
such condition is roughly matched in all cases. 



3. Results 

3.1. Comparison between the slab model and PEXRAV 

For all objects, the estimated corona temperatures (respectively optical depths) 
from PEXRAV are substantially smaller (respectively larger) than those inferred 
with an anisotropic Comptonization model (cf. Fig. la). Large differences (up 
to a factor 8) are found between the two temperature estimates. 
The reflection normalizations obtained with H94 cluster around 1, except for 
two extreme states of NGC 4151. They are, in all cases but one (ESO 141-G55), 
larger than those found using the simple cut-off power-law model (cf. Fig. lb). 
In some cases we obtain differences of factors 4-5. 

The reason of these differences is relatively simple. Within PEXRAV type 
models, the slope of the power law, determined with small errors by the LECS 
and MECS data, cannot change, by hypothesis, at higher energies. A cut-off 
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Figure 1. Slab corona temperature and reflection normalization ver- 
sus PEXRAV estimates 



around 100 keV is then required to fit the PDS data. In Comptonization mod- 
els, the LECS and MECS data determine the slope below the anisotropy break. 
Above this break the intrinsic spectrum is steeper. It can thus fit the PDS data 
without an additional steepening beyond 100 keV, allowing for a larger tem- 
perature (and consequently a smaller optical depth to keep, roughly, the same 
power law slope) and a larger value of the reflection component normalization. 



4. Discussion 

In Fig. 2a and 2b the values of r vs. kT e , obtained for the 6 sources us- 
ing the spectral models H94 and PEXRAV, are compared with the theoretical 
relations expected for a plane and hemispherical Comptonizing region in en- 
ergy balance (the temperature and optical depth then satisfy a definite and 
univocal relation which corresponds to roughly constant Compton parameters 
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with H94 (Fig. 2a) are relatively close to the 
theoretical expectations for the slab case even if they tend to fall preferentially 
above the solid line. The data therefore indicate a Comptonization parameter 
larger than for a pure slab geometry, that is a more "photon-starved" config- 
uration. The values of kT e and r obtained with PEXRAV (cf. Fig. 2b) show 
the same trend as found with H94 that is larger values of r for smaller values of 
kT e . They are also above the theoretical expectation for a slab and even above 
the theoretical expectation for a hemisphere. Therefore this set of parameters 
suggests, for each source of our sample, a configuration more "photon-starved" 
than a hemisphere. 

An important difference with respect to the slab model is that the optical depths 
derived from PEXRAV fits are generally larger than 1. The corona should then 
be optically thick reducing or canceling the effects of anisotropy. Therefore also 
the PEXRAV model has an internal consistency. However a corona with large 
optical depth may wash out discrete features from the underlying disk (e.g. lines 
and reflection itself) more than would be desirable. This problem could be alle- 
viated if the corona was " patchy" . 




Figure 2. Optical depth r versus temperature kT e from (a) H94 and 
(b) PEXRAV fits. The theoretical relations between r and kT e for 
a plane and hemispherical Comptonizing region in energy balance, in 
solid and dashed line respectively (from Stern et al., 1995), are shown 
for comparison. 



We conclude that both a hot, optically thin corona with significant anisotropic 
effects and a less hot, optically thick, patchy corona with negligible anisotropy 
are consistent with the available data. 

4.1. Correlations between physical parameters 

A correlation between the reflection normalization R and the photon index T has 
been claimed by Zdziarski et al. (1999, hereafter Z99) from the study of a large 
number of GINGA observations of Seyfert and galactic black hole objects. Yet, 
the validity of this correlation is still under debate since R and T are strongly 
correlated in the fitting procedure (cf. for example Vaughan & Edelson, 2000). 
We have plotted in Fig. 3a the reflection normalization R versus the photon in- 
dex T obtained with PEXRAV. Our data don't show a clear correlation between 
the two parameters. However our sample is biased in favor of objects with hard 
spectra, to ensure a good detection in the PDS instrument. When the 13 objects 
actually observed by BeppoSAX are taken into account a stronger correlation 
is observed (Matt, 2000), in agreement with Z99. In the case of the H94 model 
we have no simple way of characterizing the spectral shape. Also results are 
available for only 6 sources and we cannot include the other Seyfert Is observed 
by BeppoSAX. Therefore a direct comparison with the Z99 correlation is not 
possible. 

Z99 interpret the correlation between R and T in the framework of thermal re- 
processing models where R is directly proportional to the solid angle subtended 
by the cold matter surrounding the corona. In this case, the larger R and the 
larger the cooling of the corona resulting in softer spectra (larger T). 
It is interesting that the BeppoSAX data can provide the temperatures asso- 
ciated with sources with different values of R or T, which were not available 
from the GINGA data used by Z99. The result is shown in Fig. 3b where we 
have plotted Rpexrav y s- the corona temperature Ec/2. This plot suggests a 
positive correlation between the two parameters, that is, the temperature of the 
corona is larger for larger values of the reflection normalization. A plot of Eq/2 
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Figure 3. a) Reflection Rpexrav versus PEXRAV photon index T. 
b) Reflection component Rpexrav versus PEXRAV corona tempera- 
ture E c /2 from PEXRAV. The open circles are from Matt (2000) while 
the solid ones are the results of this work, c) R s i a b versus slab corona 
temperature kT e 



vs. T shows that steeper spectra correspond to higher temperatures and (as a 
corollary) to lower optical depths. 

We note that the slab model analysis suggests different trends. If we plot R vs. 
kT e for the latter model we find that sources with larger R tend to have lower 
kT e (cf. Fig. 3c). Looking also at Fig 2 it appears that this model yields a 
smaller range of optical depths for the sources in the sample, thus sources with 
flatter spectra have larger temperatures, a behavior opposite to what is found 
with PEXRAV. 

It is interesting to note that, in the interpretation of Z99, the Rpexrav~Ec/2 
correlation of Fig. 3b suggests that the corona temperature is larger for larger 
cooling. Such behavior is expected in pair dominated thermal plasma in pair 
equilibrium (i.e. where pair creation is balanced by pair annihilation). Indeed, 
in this case, an increase of the cooling corresponds to a decrease of the number 
of particles in the hard tail of the thermal particle distribution, and thus to a 
decrease of the pair production rate. The number of particles in the thermal 
bath thus decreases, i.e. r decreases. If the heating is kept constant, the avail- 
able energy is now shared among less particles, resulting in an increase of the 
temperature. 

On the contrary, the anticorrelation between R and kT e obtained with the slab 
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corona model (Fig. 3c) is naturally expected in the framework of reprocessing 
models but for low pair density corona, in complete opposition with the conclu- 
sions deduce from our PEXRAV fits. 



5. Conclusion 

The aim of this paper was to test thermal Comptonization models over the high 
signal to noise BeppoSAX observations of a sample of six Seyfert 1. We use 
two types of model: a detailed Comptonization code in slab geometry, which 
treats carefully the anisotropy effects in Compton processes, and a simple cut- 
off power law plus reflection model (PEXRAV model of XSPEC). The latter is a 
relatively good approximation to Comptonization spectra in isotropic geometry 
and/or optically thick corona. The main results of this work are the following: 

• The data are well fitted by both models and there is no statistical evidence 
for a model to be better than the other. Both models give results in agreement 
with a X-ray source geometry more "photon starved" than the slab case. 

• There are strong differences between the best fit values of the tempera- 
ture and optical depth of the corona and the reflection normalization obtained 
with the two models. We generally obtained larger corona temperature, smaller 
optical depth and larger reflection normalization with the slab geometry in com- 
parison to PEXRAV 

• The two models leads to strong different relationships between physical 
parameters. For instance, we obtain a correlation between the reflection nor- 
malization R and the corona temperature with PEXRAV, and an anticorrelation 
with the slab corona model. This has major consequences for the physical in- 
terpretation of the data. 

Forthcoming observations with CHANDRA, XMM-Newton and INTEGRAL are 
expected to bring substantial progress to discriminate between these two types 
of models for the high energy continuum of Seyfert galaxies. 
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